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Abstract 

(— I . Using the multipolar expansion of the electrostatic energy, we have characterized the asymp- 

totic interactions between an oxygen atom 0( 3 P) and an oxygen molecule 02( 3 S~), both in their 
electronic ground state. We have calculated the interaction energy induced by the permanent 
electric quadrupoles of O and O2 and the van der Waals energy. On one hand we determined the 

27 electronic potential energy surfaces including spin-orbit connected to the 0( 3 P) + 02( 3 S~) 

"X ■ 

dissociation limit of the O-O2 complex. On the other hand we computed the potential energy 
curves characterizing the interaction between 0( 3 P) and a 02( 3 S~) molecule in its lowest vibra- 
tional level and in a low rotational level. Such curves are found adiabatic to a good approximation, 
namely they are only weakly coupled to each other. These results represent a first step for modeling 
the spectroscopy of ozone bound levels close to the dissociation limit, as well as the low energy 
collisions between O and O2 thus complementing the knowledge relevant for the ozone formation 
mechanism. 
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I. INTRODUCTION 



The ozone molecule plays a crucial role in the physics and chemistry of the Earth atmo- 



sphere. However, a lot remains to 
thought to take place in two steps 



be understood, especially about its formation, which is 
l|. Firstly, an oxygen atom O and an oxygen molecule 
O2 collide to give a ro-vibrationally or electronically excited ozone complex O3. Secondly, 
this complex stabilizes by inelastic collision with one surrounding atom or molecule, which is 
the so-called deactivation process. However, this second step takes place provided that the 
excited complex O3 does not dissociate into O+O2 before colliding with the surrounding gas. 
Characterizing the competition between deactivation on the one hand, and dissociation on 
the other hand, is the key point in order to understand quantitatively the ozone formation 
in the atmosphere. 

In this respect, one of the most striking features of ozone physical chemistry is the un- 
conventional isotopic effects that influence the competition between stabilization and dis- 
sociation of O3 (2-6]. It became clear in recent years 3], that these isotopic effects were 
determined by the difference of zero-point energies AZPE of the O2 isotopologues in the 
entrance and the dissociation channels. If AZPE > 0, dissociation is energetically unfavor- 
able: the O3 has a higher lifetime, and so is more likely to give stable O3. On the contrary, 
if AZPE < 0, dissociation is energetically favorable, and it tends to dominate stabilization. 

The unconventional isotopic effects were very well understood in the beginning of the 
2000's within the framework of the statistical RKRM (Rice-Kassel-Ramsperger-Marcus) 
theory js, 9]. However, an adjustable parameter r] had to be added to that theory, in order 
to account for deviation from the energy equipartition theorem, after the formation of the O3 
molecule. Then, the need for first-principle studies of the ozone formation, especially based 
on quantum mechanics, became obvious. Since a full quantum treatment of the two-step 
process is beyond the possibilities offered by current computers, researchers are urged to 



focus on precise aspects of 
the influence of resonances 



;he process, e.g. highly-excited vibrational levels of O3 10Nl2|. 
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141 ]. or to use less demanding computational techniques, e.g. 



classical-trajectory 15[ or mixed-quantum-classical calculations 16]. 

The common features of all those studies is that they need a reliable potential energy 
surface (PES), at least for the electronic ground state X 1 A\. Since the formation of stable 
O3 involves a wide variety of geometries, from almost separated O and O2, to tightly bound 



O3, one actually needs a global PES. Up to now, all the published PESs share the same 
general features. In particular when they are cut along the minimum-energy path, they 
all show a change in character between the inner and the asymptotic regions, which is due 
to an avoided crossing with an excited electronic state, also referred to as the transition 
state. However, there is still a controversy on whether this avoided crossing induces a 



potential barrier that goes above 
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20|, or on 



or below the dissociation limit 
the contrary, a monotonic evolution of the potential energy as suggested by the most recent 
ab initio calculations 21]. 



Up to now, all the articles that aimed at describing the O— O2 asymptotic region were 
based on quantum-chemical calculations of the ground-state and possibly the lowest excited- 
states PESs of O3 
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22|. In the present paper, we propose an alternative method based 



on the multipolar expansion of the electrostatic potential energy between O and O2, both 
in their electronic ground state. This method enables us to obtain, in a single calculation, 
the 27 spin-orbit PESs connected to the dissociation threshold 0( 3 P)+02(X 3 £~). The 
calculated PESs are then functions of the electric properties, i.e. multipole moments and 
multipole polarizabilities, of the separated systems. This method is valid provided that the 
electronic clouds of O and O2 do not overlap, that is for a O-O2 distance larger than 8 Bohr. 
Such asymptotic PESs can then be used as a tool to check the quality of the ab initio PESs. 



In this article, we use two complementary approaches to characterize the O— O2 long- 
range interactions. In section [Til we consider the atom and the diatom at fixed geometries, 
that is at given inter-particle distances and bending angles. The obtained PESs can be 
directly connected to ab initio PESs. In section IHIt we include in our model the vibration 
and rotation of O2 which yields one-dimensional potential energy curves (PECs) depending 
on the O-O2 distance. Such curves are relevant for the low energy dynamics of the O- 
O2 complex, in particular to which extent the rotation of O2 is hindered by the presence 
of O. These curves are also useful for modeling the vibrational levels of ozone close to its 



dissociation limit 1CH12| . We discuss in section IIVI the way to connect both approaches 
which is done for the first time, to the best of our knowledge. Section [V] contains our 
conclusions and prospects. 



II. ASYMPTOTIC ELECTRONIC POTENTIAL ENERGY SURFACES BETWEEN 
O AND 2 



As the two ground state particles O and 2 are far away from each other, i.e. their 
electronic clouds do not overlap, we use the Jacobi coordinates to describe the PESs. The 
body-fixed frame associated to the O-O2 complex has its z axis connecting O to the center 
of mass of O2. The x axis is perpendicular to the z axis and is located in the plane of the 
three atoms; the y axis is perpendicular to this plane. We also introduce the coordinate 
system Xdy d Zd linked to the O2 diatom. The axes Zd and Xd are related to z and x by a 
rotation of angle 9 about the yd = y axis. We denote as r the interatomic distance in the 
O2 molecule, and R the distance between the O atom and the center of mass of O2. The 
formalism presented in this section allows for calculating three-dimensional PESs depending 
on R, r and 9. In the following, we set r = r e = 2.282ao («o = 0.0529177 nm is the Bohr 
radius and the atomic unit for distances) at the equilibrium distance of O2 without loss of 
generality for our purpose. 



A. Model 



The calculations below are based on the formalism described in Ref. 23|, and used in 



Rcf. 



24j. The quantities referring to O and O2 are respectively characterized by the sub- 



scripts "a" after "atom" and "d" after "diatom". The starting point of our model is the 
multipole expansion in inverse powers of R of the electrostatic energy between two charge 



distributions "a" and u d" (see 



23|, Eqs. (2) and (3)) 



(i) 



m=— £< 



with = mm(£ a ,£ d ) and [25 1 



Qmk^ai ^d) 



-1)^(4 + 4*)! 



V(L + m)!(4 - m)\{£ d + m)\{t A - m)\ 



(2) 



In Eq. ([T]), Q™ a (Q t ™) is the i a - (£4-) rank and m-(-m-) component multipole-moment 
operator of the atom (diatom) expressed in the xyz coordinate system. The multipole 



moments of the diatom is conveniently expressed as 

Q1T= E d -mm d (m7 d d , (3) 

™, d =-e d 

where q^ d is the multipole moment operator of the diatom expressed in the XdUd^d coordinate 
system, and dt d mmd {0) is the reduced Wigner function, which characterizes the rotation from 
the XdUdZd to the xyz frame. 

The multipolar expansion ([T]) is valid provided that the electronic wave functions of O 
and O2 do not overlap. This condition is fulfilled for distances R larger than the so-called 



LeRoy radius 26| Rlr = 2{-y/ (r 2 )o + a/ (r 2 )o 2 } ~ & a o, where the mean squared radius of 
the electronic wave function in O is (r 2 ) Q = 2.001a . The quantity (r 2 )o 2 ~ r e /2 + (r) Q 
is roughly evaluated from the equilibrium distance r e = 2.282ao of O2 and the mean radius 
of the electronic wave function in O (r)o = 1.24ao- 

In this work, we calculate the two leading terms of the R~ n expansion for R > Rlr- the 
first-order term reflecting the interaction between the permanent quadrupole moments of O 
and 2 (£ a = = 2) scaling as R~ 5 , and the second-order (Van der Waals) term related to 
the interaction between the induced dipole moments (l a — Id — 1) scaling as R~ 6 . While 
generally tedious to obtain from ab initio calculations higher-order contributions may also 
be significant around Rlr- For instance, assuming that the maximal values of the Cq and Cg 
coefficients for the O-O2 interaction are half of those for O2-O2 {27]]) we obtain Cq ~ 31 a.u. 
(close to the values obtained in the present work, see next section) and C§ ~ 1291 a.u.. At 
R = 8ao the Cg term would then represent at most 64% of the Cq one and 41% at R = 10ao- 
Therefore our work provides the essentials of the long-range interaction between O and O2. 
In order to match the present asymptotic expansions to the ab initio calculations, the safest 
way is to use the C5 and Cq values determined below to fit the long-range part of the ab 
initio PESs beyond Rlr and to extract the related higher-order terms. Another possibility 
would be to directly compute the next term Cg, but this beyond the scope of the present 
paper. 

1. Zeroth-order energies and state vectors 

The oxygen atom is in an arbitrary fine-structure level \J a M a ) of its ground state 3 Pj a M a , 
where J a is the total angular momentum with projection M a on the z axis, resulting from 



the sum of the projections M^ a and M$ a of the orbital L a and spin S a angular momenta 
of the atom, respectively. Since the spin-orbit constant Ao of oxygen is negative (Ao = 
—79.1 cm -1 ), the 3 P2 fine-structure level is the lowest in energy. 

The ground electronic state 3 £~ of O2 has an orbital angular momentum projection 
Ad = on the Zd axis and a spin Sd = 1 with projection on z d . The fine structure in the 
O2 energy spectrum induced by the spin-spin interaction writes 



AdSd^d 



H. 



ss 



AdSd^d 



A 
2 
A(r 



2>S 2 Z - S 2 



(3Z d -S d (S d + l)) 



(4) 



where A(r) is the (r-dependent) spin-spin constant, with A(r e ) = 1.980 cm -1 28|. The 
interaction between the electric multipole moments of O and O2 only depends on the spatial 
coordinates of the electrons and the nuclei so that we will ignore the O2 fine structure in the 
following. Assuming an energy origin at the 3 P2 oxygen level and at the O2 ground level, 
the zeroth-order energy is 



En 



A 



o 



{Ja (Ja + 1) " 6) 



(5) 



corresponding to the unperturbed state vectors |^ ) = \JaM a ) \A d T> d ). 



2. First-order quadrupole-quadrupole interaction 



The first-order quadrupole-quadrupole interaction is obtained by setting £ a 



2 and 



n = 5 in Eq. (3) of Ref. [23|. Since L a = S a = 1, we can rewrite Eqs. (4)- (5) of Ref. [23| as 



j^JaM a J' a M' a k d T. d k' d T.' d ^ ^ ^ 

1 A \^ r J a M a 

^E 6 ^K 1^ C lM La lM Sa 
2^ ^lM'lM' Sn °M Sa M' Sa 



M' T ML 



x njj V 5 f^ (r)d 2 _ m (8) 



(6) 



■1C-, 



where refers to the components of the 2-rank tensor in the diatom frame, C c aabf3 — 
(aab/3\cj) is the compact expression of Ref. 29j for the Clebsch-Gordan coefficients, and 
where 

ii,jj 'V^ md (r) = gULJ d )(lM L JQ^\lM' La )(Ad\qT d (r)\A' d ) ■ (7) 
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The superscripts and j,f designate in a compact way the quantum numbers of and 
2 , respectively, i.e. i = {L a = l,M La }, i' = {L' a = 1,M' L J, j = {A d = 0, E d } and 
j' = {A' d = 0, £ d = Erf}. As the 2 electronic ground state is of £ symmetry the only non- 
zero component of the quadrupole operator is q(r) = (Ad = 0| q™ d ~° |A^ = 0). In Eq. ([7]) 
we only keep the quantum numbers which are not fixed: u ' JJ '^ 2 mm and E^ IaJ °- MakdT ' dk <i T 'd 
can be simplified to MLaM ' L -V 5 ^ md and E^ aJ ' aM ' a respectively 

Using the Wigner-Eckart theorem, we can connect all matrix elements of O quadrupole 
moment to a single one, say that for which Mi a = M' L — m — 0, 



<im^q-|imO = (-i) M£ "Vt( 1 2 1 

V 2 ^ -M La m M' La 



Q, (8) 



with Q = (10|Q 2 |10), and where the symbol (...) is a Wigner 3-j symbol, which imposes 
M L = M'r + m. 



3. Second- order dipole-dipole interaction 



The second-order interaction scales as R~ 6 and results from the dispersion term due to the 
induced dip ole- induced dipole interaction. It is calculated by setting £ a = i' a = t d = £' d = 1 
and n = 6 in Eqs. (6)-(10) of Ref. 23|]. In this case, the coupled polarizabilities (Eq. (7) 
of Ref. 23|) are tensors which can be of rank C a = 0, 2 and £d = 0, 2, for O and 2 
respectively. Then, the dispersion energy (see Eq. (8) of Ref. 23] ) can be written as a sum 
of an isotropic (C d = 0, i.e. ^-independent) and an anisotropic contribution (£ d = 2, i.e. 
dependent) 

M La M Sa 

E r J 'a M a X 
°1MC 1M' °M Sa i 



W JaM a J'aM a( p 

-^disp I'l^l 



X 



, M> 



Ml ML 



x 



+ £ 



2 



M-- 



C< M d =-2 



1 j T/ 

La V 



62MM d 



(r)d 2 _ MM (9) (9) 



with £ < = mm(£ a , £ d — 2) = £ a - Note that in Eq. (jUJ), the labels for E disp and V dlsp have 
been simplified in the same way as for the first-order term (see text after Eq. (j7J)). 

The dispersion coefficients MhaM La V^ MMd (given by Eq. (9) of Ref. 23]) depend on dy- 
namical polarizabilities at imaginary frequencies. They are conveniently expressed in terms 
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of coupled polarizabilities for O and O2 respectively, following the definitions introduced in 
Ref. |3o| 

M ^ M '^a {11)CaM (icu) = ^(lmlm'|£ a M) M ^ M ^a lmW (zo;), (10) 

mm' 

ot(u)c d M d (wj] r ) = ^2(lqlq'\C d M d )a lqlq/ (iu; r) (11) 

qq> 

The uncoupled dynamic polarizabilities Mia ctioio = Mha a zz for the M La sublevels of 0( 3 P) 
assuming an electric field polarized in the z direction 



iave been calculated by one of us using [N,N-1] Pade approximants ( La ak) and ( La bk) 



31] reported in the supplementary material [32| for convenience, and published elsewhere 
33|. For 0( 3 P) we obtain 



u La M' ( , lA M ^=°a zz (iu;) + 2 M ^= ±1 a zz (zu J ) 

L *a {ll)m {iu) = -j= (13) 



ML ° M '^a {n)2M (icu) = (-l) 1 "^- VE 



1 2 1 

-M La M M' La 
x ( M ^ =±1 a zz (iuj) - J **- =0 a„(iw)) . (14) 

Similarly the coupled dynamical polarizabilities of O2 ct(n)c d M d are related to the (r- 
dependent) uncoupled ones a ww (iu;r) = a"(iuj;r) (the parallel component along z d ) and 
Qj-in-i (icu; r) = —a L {iuj;r) (the perpendicular component with respect to Zd) according to 



23] 



a^(iu;r) + 2a ± (iu;r) 
«(ii)oo(2w;r) = (15) 

a(ii)aMi(iw; r) = ^|<W (a^iw; r) - a L (iu\ r)) . (16) 

Note that just like (see Eqs. (jSJ) and (JTj)), A^d is zero because O2 is in a £ electronic 
state. 



B. Asymptotic PESs: results and discussions 

We have calculated the permanent quadrupole moment Q of 0( 3 P) (see Eq. (JE}), with 
the CASSCF method, in a full-valence active space including 6 electrons and 5 orbitals 
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TABLE I. Quadrupole moment Q and static dipole polarizabilities MLa a zz (uo = 0) (in a.u.) for 
0( 3 P). All values are expressed in atomic units. 

Ref. Q ML ° =0 a zz M ^ =1 a zz 

This work" -0.95 5.64 4.83 
Ref.[34f -0.95 

Ref.[35j c - 5.86 4.94 

Ref.[36j d -1.02 5.91 4.89 

Ref.[36j e -1.04 6.08 4.99 



a Full valence CASSCF(6e,5o) with aug-cc-pVQZ basis set 
& CCSD(T) with aug-cc-pV5Z basis set 
C CCSD(T) with quadruple-C GTO/CGTO basis set 
d CASSCF and e CASPT2 both with triple-C GTO/CGTO basis set 

and an "aug-cc-pVQZ" basis set. Our value Q = —0.95 a.u., is in very good agreement 



with the value obtained with the CCSD(T) using quadruple-^ (or higher) basis set 34J of 
similar quality to the one presently used (Table H)). This indicates that the contribution 
of dynamical electron correlation effects which are missing in our CASSCF treatment can 
safely be neglected, provided that sufficiently large basis sets and active space are used. 
Only static dipole polarizabilities (u = 0) have been reported up to now in the literature, 
and Table [J shows that our value obtained from Eq. (|T2|) is in satisfactory agreement with 
other published values. 

The available values of the permanent quadrupole moment and static dipole polariz- 
abilities of 02(X 3 S~) are shown in Table HT1 In the present work, the quadrupole moment 
q(r) = (Ad = 0| g° I Ad = 0) is taken from Ref. [s^], where it is calculated at the equilibrium 
distance r e = 2.282 a.u. and for the lowest vibrational level v a = using an harmonic- 
oscillator approximation. The dynamic dipole polarizabilities are taken from the semi- 



empirical dipole-oscillator-strength distribution (DOSD) values of Ref. 38|. All quantities 



27] 



regarding O2 are in good agreement with the best recommended ab initio values of Ref. 
as well as with experimental values. 

Tables UTTl and HVl present the long-range multipolar coefficients at the first (electrostatic) 



TABLE II. Quadrupole moment q, and parallel a^(ui = 0) and perpendicular cr^w = 0) static 
dipole polarizabilities (in a.u.) for 02(X 3 £~). All values calculated ab initio are given at the 
equilibrium distance r e = 2.282 a.u. and are expressed in atomic units. The value in parenthesis 
is calculated for for the lowest vibrational level following Ref. 371 ] . 



Method 



Qjll 



a 



Ref. [37] a 
Ref. [38] b 
Ref. [27] c 
Experimental 



-0.2530 (-0.2273) 

-0.2251 
-0.3±0.1 d 
-0.25^ 



15.29 
15.367 
15.7±0.3 e 
15.379 



8.24 
8.228 
8.4±0.3 e 
8.229 



a CBS-CASSCF+1+2 
6 semi-empirical DOSD values 
^Recommended values: CAS(12e-14o)-ACPF calc. with aug-cc-pV5Z basis set 



39] 



^optical birefringence 
e vibration rotation Raman spectroscopy 
' pressure-induced far-infrared spectrum 

9 depolarization ratios 42] 



4l| 



TABLE III. Long-range non-zero Ml ° Mi « V^^ electrostatic coefficients (in atomic units, see 
Eq. d2J) for the 0( 3 P) + 2 (X i T l ~ ) interaction at r = r e = 2.282 a.u. and for the lowest vi- 
brational level Vd = using the data of Table [III 



M La Ml m Vt c (r e ) Vt c (v d = 0) 

±1 ±1 -0.721 -0.647 

±1 ±1 0.832 0.748 

±1 Tl ±2 -0.294 -0.264 

±1 =Fl -0.832 -0.748 

1.442 1.295 
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TABLE IV. Long-range non zero Ml ° Ml « VqCjMO dispersion coefficients (in atomic units, see 
Eq. ©) for the 0( 3 P) + 2 (X 3 S-) interaction at r = r e = 2.282 a.u.. 



M La Ml M U V^T 

±1 ±1 -30.24 

±1 ±1 2 -3.665 

±1 ±1 2 -0.253 

±1 Tl ±2 2 0.179 

±1 Tl 2 0.253 

-31.511 

2 -4.322 



and second (dispersion) orders of the perturbation theory for the matrix elements of O + 
2 . The ML " M ' L "Vg% and MLaM ' L *V^ M0 coefficients are given for each (M LaJ M' La ) pair 
and for various values of m, and M. at r = r e = 2.282 a.u., and also when the relevant 
quantities are averaged over the Vd = level of O2. for the former. As already mentioned 
previously, the maximum value of the dispersion coefficients are about half of those for the 
O2+O2 long-range interaction of Ref. j^] which represents a good test of the consistency of 
the calculations. 

Asymptotic PESs are obtained after diagonalizing the total interaction potential matrix 
with elements 

E t tf M (r, R, 6) = 5 JaJ ,5 MaM ,E + E^ K (r, R, 6) + i^«(r, R, 9), (17) 

in the subspace spanned by J a and M a . Figure Q] displays one-dimensional cuts of these PESs 
at r = r e = 2.282 a.u. either for a fixed bending angle (8 = 0, FigJUa)) or at a given O-O2 
distance (R = 10 a.u., Figfjjb)). In Fig.[T^b) we note that the first two states degenerate into 
a IT state for collinear arrangements (#=0 and 180°). Dispersion contributions are noticeable 
in asymptotic ozone and change slightly the anisotropy of the ground-state potential which, 
however, remains almost isotropic. The second A' state becomes attractive after inclusion 
of dispersion energies. 

Inclusion of spin-orbit splitting lifts the degeneracy of the states and twenty-seven states 
arise from the interaction of 0( 3 Pj a ) with 02(X 3 S~), which reduce to nine states if we neglect 
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FIG. 1. Asymptotic PESs for the three electronic states 1A' (solid line), 1A" (dashed-dotted line) 
and 2A' (dashed line) correlated to 0( 3 P) + 2 (X 3 S~) for r = r e = 2.282 a.u.. (a) at 0=0°; (b) at 
R = 10 a.u.. In (b), the curves with filled circles refer to the PES when only quadrupole-quadrupole 
terms are included. 



the O2 fine-structure (Fig. Eta))- The angular dependence of the five states correlated to 
the 0( 3 P2)+02(X 3 S~) limit are displayed for R = 10 a.u. on Fig. EJb) , where the influence 
of the spin-orbit interaction of 0( 3 P2) on the anisotropy of the PES is clearly visible. 
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FIG. 2. (a) Asymptotic PESs for the nine spin- orbit states O( 3 P ,i, 2 ) + 2 (X 3 £~) at r = r e = 
2.282 a.u. and 9=0°. (b) Asymptotic PESs for the five lowest spin-orbit states correlated to 0( 3 P2) 
+ 0(2X 3 S~) at r = r e = 2.282 a.u. and R = 10 a.u.. In (b) the dotted line refer to the ground 
state PES without fine-structure. 

III. ASYMPTOTIC POTENTIAL ENERGY CURVES BETWEEN O AND 2 IN 
A GIVEN RO VIBRATIONAL LEVEL 



A. Model 



The calculations presented in this section are based on our previous work on Cs(6p 2 P)- 



Cs 2 ( 1 E+ v d 



0,N d 



25 



43 



45| performed in the context of ultracold gases. The main 



differences are that the present work deals (i) with a ground-state atom, and (ii) with a 
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triplet molecule (Sd — 1). 

When an oxygen atom approaches a rovibrating O2 molecule from large distances, the 
rotational levels of O2 are coupled by the electric field induced by O. We focus on the deriva- 
tion of one- dimensional long-range potential energy curves (PECs) for the O-O2 complex 
depending on R, thus describing the interaction in the frame xyz linked to the complex. In 
other words we leave out here the mutual rotation of O and O2. 

The quantities regarding the oxygen atom are expressed here in the fine-structure basis 
\J a M a ) and are identical to those of Section [Til To stress the differences with the fixed- 
geometry problem for the diatom, all the quantities X that are averaged over O2 vibrational 
wave functions are overlined (X). Note that the effects of centrifugal distortion on those 
wave functions will be ignored. 



1. Zeroth- order energies and state vectors 

We consider that the O2 molecule lies in its ground vibrational level v d = 0. In addition 
to the electronic quantum numbers = 0, Sd = 1 and of O2, we introduce the quantum 
numbers for the (electronic+nuclear) orbital angular momentum Nd, for the orbital+spin 
angular momentum J d (\N d — 1\ < J d < A^ + l), and the corresponding projections M d , M Nd 
(and Ms d for Sd) on the z axis. Here we limit our study to the 16 C>2 isotopologue allowing 



for odd values of N d only 46] . All the matrix elements will be given in the fine-structure 



basis \(N d Sd)JdM d ) = \N d J d M d ), connected to the \N d M Nd S d M Sd ) basis by 



\N d J d M d )= C N d M N ,s d M Sd \N d M Nd S d M Sd ) 



(18) 



M Nd M Sd 



where C J J K Ii 



N d M Nd S d M Sd 



{N d MN d SdMs d \JdMd) denotes a Clebsch-Gordan coefficient. 



In its ground electronic state, O2 belongs to Hund's case b, i.e. Nd is considered as a good 
quantum number. The rotational spectrum is dominated by the free-rotator contribution 



B Vd N 2 d , with B 



v d =0 



B n = 1.438 cm" 1 



46j . The matrix element of the spin-spin interaction 
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(Eq. are expressed in the | J d M d ) basis (see Sec. Ill of the supplementary material [32|) 



E ss = (J d M d \H ss \J d M d ) 



3 



(2N d + 1) A 



o 



N d 1 J d \ N d 2 N d 
1^ 2 



where the spin-spin coupling constant Ao = A„ d= o = 1.983 cm -1 is taken from Ref. 28|. The 
spin-spin interaction results into the splitting of the fine-structure rotational levels J d inside 
a given manifold N d . The spin-rotation interaction Hsr = 7^ Vd= o^d ■ is characterized by a 
coupling constant Jt Vd= o = 8.43 x 10 -3 cm -1 28| much smaller than Ao, and will be neglected 



in what follows. 

The zeroth-order energy Eq is the sum of the atomic spin-orbit, the molecular rigid- 
rotator and the spin-spin interactions 

E = \ (Ja (Ja + 1) - 6) + B (N d (N d + 1) - 2) + E ss , (20) 

corresponding to unperturbed state vectors j^o) = \J a M a ) \N d J d M d ). The origin of energies 
E = is fixed to the 0{ 3 P2)+02(v d = 0, N d = 1) dissociation limit. 



2. First-order quadrupole-quadrupole interaction 



The quadrupole moment matrix elements of the vibrating and rotating O2 molecule is 
obtained by starting from Eq. (12) of Ref. 44J written in the \N d M^ d S d Ms d ) basis, and by 
applying the transformation to the fine-structure basis \NjJdMd) (Eq.( ll8j) ) 



(v d = 0, N d J d M d \Qz m \v d = 0, N' d J' d M' d ) 



4) i+wj-*s- / (2 j d+ i) ( 2 J' d +l) 



xyJ(2N d + l)(2N' d +l) 



N d 2 N' d 




K 1 J' d 



2 J' d 



(21) 



Jd 

J d 2 N d J y -M d -m M' d 
where we use the expansion of the product of three Clebsh-Gordan coefficients (see Eq. (2) 



of the supplementary material 



32|), and where 

+00 



% — Qv d =o 



dr (ijj Vd =o(r)) q{r) 



(22) 
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is the quadrupole operator averaged over the vibrational wave function tf) Vd= o(r). Its value 
q = —0.2273 a.u. is taken from Ref. [3?]]. The properties of the 3-j and 6-j symbols impose 
that N' d = N d , N d ±2 and J' d = J d , J d ± 1, J d ± 2. 

The independent quadrupole-quadrupole matrix element reads 

-^J a M a J' a M' a N d J d M d N' d J' d M' , 
1 2 

= & E 9m(2,2)(J a M a \Q^\J a M' a ) 

m=-2 

x (v d = 0, N d J d M d \Q~ 2 m \ Vd = 0, N' d J' d M' d ) , (23) 

where g m (2,2) is given by Eq. (T5]). Note that the atomic quadrupole moment, expressed in 
the fine-structure basis, is related to the one in the \M La M Sa ) basis (see Eq. ([H])), 

(j a M a \Q™\ j' a M' a ) = y] y] c ^t a iM Sa c iM' L \M' Sa 

M La M Sa M' La M' Sa 

xS MsaMk {lM La \Qf\lM' La ). (24) 



5. Second-order dipole-dipole interaction 



In order to calculate the dispersion term in R~ 6 , we use the approach of Refs. 43|, l44j . 
adapted with the notations of the present paper. In the atomic and molecular fine-structure 
bases, the matrix elements associated with the second-order dipole-dipole interaction reads 



-=J a M a J' a M a N d J d M d N' d J d M d , 
-^disp I-") 

! 1 



2ttR 6 



m,m'=— 1 



x / du JM ^a lmlm \iu) 



x^ J -^-«a 1 _ ml _ m ,(i w ) J 



(25) 



where JaMaJ ' aM ' a ai m i m >(iui) and NdJdMdN 'd- J 'd M 'dai- m i- m '(iu}) are the uncoupled dipole polar- 
izabilities at imaginary frequencies, for the atom and the molecule respectively. Note that 
these polarizabilities are related to the ones defined in Ref. 43] according to xx ' a\ m \ m i{ioj) = 



[(— l) m } xx 'a^ xA ^,(iu), where x stands for the set of atomic or molecular quantum numbers 



(see Ref. 



Eqs. (13) and (14)). The polarizabilities of O expressed in the fine-structure 



1(3 



basis reads 



lml m '(«w) =^ C 

M La M Sa M' La M' Sa 
X 6 



JaM a r J' a M' a 

lM La lM Sa ^lM' La lM' Sa 



>M Sa M' s 



Sa 



ML * M '^a lmlm ,{iu) . 



(26) 



where MLaML aai m i m >(iuj) are the polarizabilities in the \M La Ms a ) basis (see Eqs. (TT3~j) and 

(HI)- 



material 



For Q 2 , the polarizabilities JiMdJ ' dM ' d oi\ m \ m i are calculated in Sec. IV of the supplementary 



32], starting from our previous work 43j, |44j. We obtain finally 



N ^ M ^<a lmlm ,{iu) = {-l) m 5 JdJ>d 5 MdM>d 5, 



m,—m' " 



al(iu) + 2«q (iu) 



+ J H A /(2iV d + 1) + 1) (2 J, + 1) + 1) 



^ 1 ^ 

2 iV d 



^ 2 




1 1 



J,/ 



4 



M' d -M d m m' ) \ -M d M d - M' d M' d 



(27) 



where a{j = ol" d=0 and o7q = ®t d =o are the parallel and perpendicular polarizabilities of 2 
in its vibrational ground level. Equation ( 12?]) is a sum of two contributions: The first term, 
proportional to the so-called isotropic polarizability (ajj + 2o7q )/3, is diagonal; the second 
term, which is proportional to the so-called anisotropic polarizability ajj — o7j}, couples the 
different angular-momentum projections M d . The properties of the 3-j and 6-j symbols 
impose that N' d = N d , N d ±2 and J' d = J d , J d ±l,J d ± 2. 



B. Asymptotic potential energy curves 

The asymptotic PECs are obtained after diagonalizing the potential energy matrix with 
elements (where we omitted the label v d = for simplicity) 

^J a M a J a M a N d J d M d N' d J' d M' 

-^tot [ K J - d JaM a O,J' a M' a Oj d M d Oj' d M' d ^0 

^J a M a J' a M' a N d J d M d N' d J' d M' 

-^J a M a J' a M' a N d J d M d N' d J' d M' 
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Atom-diatom distance R(a.u.) 



FIG. 3. Long-range potential energy curves for the interaction between 16 0( 3 P2,i,o) and 
16 02(^ 3 S~, Vd = 0,Nd = 1,3,5), as functions of their mutual separation R, for M = and 
the (+) reflection symmetry. 



for different values of R, and within the subspace determined by each value of the total 

B,^]- For M = 0, the eigenvectors 



43 



angular momentum M = M a + Md on the z axis [25 , 
of ( )28|) are also characterized by a given reflection symmetry through any plane containing 
the z axis. Following the expression of the reflection operator a xz for two atoms (see for 
instance Eq. (3.11) of Ref. [48]), we assign the ± symmetry to eigenvectors corresponding to 
the linear combinations \M a , -M a ) ± (-\) L «+ J *+ N d+Jd\ _ M a M a ). On Figs. |3H6l we display 
long-range potential curves belonging to the M = + symmetry. The curves belonging to 
other symmetries have the same appearance, except that some asymptotic channels may not 
be allowed for a given value of M, e.g. J a = Jd = for M — 1. All the curves are strongly 
attractive, due to the dominance of the attractive van der Waals term in the range of O-O2 
distances that we consider here (see Tables HTT1 and HVl for an illustration of this point). 

As the spin-orbit splitting of O is much larger than the rotational splitting of O2 (at least 
for the Nd lowest levels), we see in Fig.|3]that the corresponding fine-structure PEC manifolds 
are almost decoupled from each other. This situation corresponds to Hund's-like case " 1C" 
defined in [4^], also valid in the Cs*+Cs2 system 45]. For the sake of clarity, we have 



restricted our plot to the first three (odd-parity) rotational levels of O2 (Nd = 1, 3, 5), even 
if in principle all the rotational levels up to infinity should be included in the calculation (see 
a discussion about this point in 44j). However, even at the \ow-R limit of our calculations, 
the different Nd levels remain almost uncoupled. As an illustration, at R = 8 a.u., the PECs 
of Figs. [3H5] conserve at least 95% of their asymptotic N d character. In other words, in the 
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8 10 12 14 16 18 20 

Atom-diatom distance Ft (a.u.) 



FIG. 4. Same as Fig. O but focused on the 3 i-*2 fine-structure level of 16 0. 




Atom-diatom distance Ft (a.u.) 

FIG. 5. Same as Fig. El but focused on the 3 i"2 fine-structure level of O, and the Nd = 1 rotational 
level of 16 2 . 

range of O-O2 distances that we consider, the oxygen molecule rotates almost freely. 

The PECs obtained for other O2 isotopologues, not shown here, have similar features 
as Figs. [3H5j But caution has to be taken on which Nd values are allowed. For 18 O2, the 
allowed N^ values are the odd ones, whereas for 17 02, the allowed values are the even/odd 
ones, if the total nuclear spin is odd/even. For mixed isotopologues, that is 16 17 0, 16 18 
and 17 18 0, all values of N d are possible. But since the Hamiltonian of Eq. (128]) conserves 
the parity of Nd, even and odd values can still be separated from each other. 

In order to get an insight into the collisional dynamics between O and O2, we have added 
to our curves a centrifugal term L(L + l)/2fiR 2 , where L is the partial wave of the collision 
that we took equal to 10, and \x the reduced mass of the O-O2 system (Fig. Ej). This term 
competes with the attractive van der Waals interaction, and creates a centrifugal barrier of 
about 2 cm -1 , or about 3 Kelvin. This temperature being much lower than even the lowest 
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10 



Atom-diatom distance Ft (a.u.) 



100 



FIG. 6. The long-range potential curves of Fig. [5] (solid lines) dressed with a centrifugal term 
corresponding to a collisional partial wave L = 10 for the 16 0"2- 16 system. Note that the R axis 
is logarithmic. 

temperatures of ozone samples, we can conclude that many partial waves should be included 
in a full quantum treatment of the O-O2 collision. 

IV. DISCUSSION: HOW DO THE TWO APPROACHES CONNECT TO EACH 
OTHER? 

In the previous sections, we have calculated the matrix elements of the potential en- 
ergy operator of the long-range interaction between ground state O and O2, yielding its 
eigenvalues after diagonalization in the appropriate configuration subspace: either for fixed 
geometries of the O-O2 complex leading to (R,9,r = r e )-dependent electronic PESs (Sec- 
tion Hi]), or for a given rovibrational state of O2 leading to independent PECs (Section HID . 
It is clear that those two adiabatic potential energy sets are not directly related to each 
other; instead, the matrix elements of the latter approach are obtained by averaging those 
of the former approach on the r and 9 coordinates of O2, as briefly addressed in Ref. {47) ]. 
This is easily demonstrated on the matrix elements of the quadrupole-quadrupole interac- 
tion. Equation ( 1211) shows that the ro-vibrationally averaged matrix element of the diatom 
quadrupole moment reads 



( 



v d N d J d M d |g 2 " m | v d N' d J' d M' d ) = q. 




'7T 



d9 sin 9i/} Nd j dMd WtmoWN'^M'j {9) 



(29) 
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where 

^N dJdM M = \Fy^ E (N d M Nd S d M Sd \ J d M d )dX M (30) 



M Nd M Ns 



is the angular part of the wave function of the fine-structure rotational level. We thus have 

Jo Jo 

The situation for the second-order dipole-dipole interaction is slightly more subtle. In 
Sec. V of the supplementary material [32j, we demonstrate that as far as the rotation of O2 
is concerned, the passage from one approach to the other is made by averaging -Edisp over 
the rotational wave function i/jN d j d M d {@), which can be schematically summarized as follows 

d c J M0 (6) [in E disp ]^ / ^sin^iv d j d M d (^)^jwo(^)^j^(^)[inS disp ], (32) 
Jo 

where E disp = E^ M and E disp = W are given by Eqs. ® and 

(|25|) . respectively. 

The drawback which can occur when taking in account the vibration of O2 is illustrated 
by considering the static parallel electronic polarizability a" (a; = 0; r) 

a "^0;r) = 2^ (33) 

where Vx{r) and V e (r) are the potential energy curves of the 3 S~ ground state and of 
the excited electronic states, respectively, and Hxeij) the (r-dependent) transition dipole 
moment between states X and e. Eq. ( 133]) is a strict application of the Born-Oppenheimer 
approximation, i.e. the molecular response to an electric field is only due to the electrons. 
The most straightforward way to take into account the vibration of O2 would be to average 
( )33l) over the vibrational wave function, 

(Ar% d = 2^| + ^^Jr)^|^^(r). (34) 

However, one sees that Eq. (1341) differs from the polarizability a|J associated with the vi- 
brational level v d , which reads 

e,v^X,v d ^ ^ 

where E ev i and Ex Vd are the electronic and vibrational energy levels of the field- free molecule. 
Note that v' d can also stand for a continuum state, in which case the discrete sum in Eq. ( 135]) 
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becomes an integral. As pointed out in Ref. |49|, Eqs.( l34p and (|35|) are identical if one 
neglects the diatom zero-point energy, i.e. Ex Vd = Ex(r = r e ) and E ev r = E e (r = r e ), 
where r e is the O2 equilibrium distance in its ground state. 

To check the validity of this approximation is appropriate for O2, we have considered the 
two lowest excited B 3 Y,~ and E 3 H~ electronic states in the sums above (Fig. [7]). The X 
and B PECs, as well as the transition dipole moment fixB(f) are taken from Refs. 5(], 51 1. 



The PEC for the pure E Rydberg state is taken from Ref. 52J, and we have extended it 
for 2.8ao < r < 3ao by shifting down the experimental potential curve of Ot(X 2 Il g ) 53]. 



The transition dipole moment between X and E has been extracted from Ref. [52J. As 
described in Ref. [54j . the B and E states exhibit an avoided crossing around 2.3ao close 
to the X equilibrium distance r e , due to a valence- Rydberg exchange of character. The B 
state is steep and purely repulsive around r e , while the minimum of the E state is located 
at r = 2.25ao; namely close to r e . 

The numerical values of Eqs. (I33H35P obtained within this three-state model are presented 
in Table IVl The electronic polarizability of the X state a"(r e ) =5.75 a.u. amounts for 40% of 
the total parallel polarizability 38j. The three determinations are almost identical, despite 
the particular shape of the PECs. The partial contribution of each transition, (B 4— X) 
and (E <— X) show a discrepancy between a"(r) on the one hand, and (a"(r)) and a{j 
on the other hand, which is due to the sudden changes in potential curves and transition 
dipole moment within the spatial extension of the X, = level. The fact that these 
differences compensate each other is due to the exchange in character between the two 
excited electronic states. In contrast, inside each electronic transition, (a"(r))„ d and a| 
are always nearly equal. On a physical point of view, it confirms that the response to an 
external electric field is, to a very large extent, due to the electrons, and not to the nuclei, 
which justifies the use of the Born-Oppenheimer approximation. 

Our result tends to prove that the vibrational part of the energy (see Eq. (|35|) ) can be 
neglected, whether the excited state PEC has a potential well or is purely repulsive. In the 
first situation, which corresponds to the E 4— X transition, the ground vibrational level 
Vd = has significant Franck-Condon factor (higher than 10~ 4 ) with the E vibrational lev- 
els ranging from v = to 4. The vibrational energies (with respect to the corresponding 
minimum electronic energies) associated with X, va = and E, v = 4 are 791 cm -1 and 
10300 cm -1 respectively. So they are very small compared to the difference in electronic en- 
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FIG. 7. Potential energy curves of O2: (a) for the B 3 T,~ state (dashed lines) and the -E 3 S~ state 
(dotted lines), (b) for the X 3 Ti~ state, (c) transition dipole moments: B X (dashed lines), 
E X (dotted lines). Between the two vertical lines, the wave- function of the X, = level is 
more than half its maximum value. 



TABLE V. Static parallel polarizability of O2 calculated with three electronic states (see text): 
a'l(r e ) is purely electronic polarizability, evaluated at the O2 equilibrium distance (see Eq. (|33|) ). 
(a"(r)) t , d= o is the r-dependent electronic polarizability averaged over the ground-level vibrational 
wave function (see Eq. (fM|) ). and a?L = is purely vibrational (see Eq. (f35j) ). We present the 
contribution of individual electronic transitions, B <— X and E <— X, to the polarizability. 

Transition a^(r e ) (a 11 (r)) Vd=0 ®l d =o 

B^X 5.47 5.16 5.15 

E <- X 0.29 0.63 0.63 

Total 5.75 5.79 5.78 
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ergy, which equals to 80000 cm" 1 at r = r e . The ground vibrational level Vd = significantly 
overlaps with continuum states of the purely repulsive B state whose classical turning point 
is located in the Franck-Condon region. This turning point region is precisely where the 
vibrational, i.e. kinetic part of the energy is small. In the present case, we only have the O2 
polarizabilities at the equilibrium distance r e . In consequence, we assumed that the rela- 
tionship between the fixed-geometries and the ro-vibrationally averaged approach, given for 
-Eeiec in Eq. fl3T|) . can be extended to -Edisp, and hence for E tot , within a good approximation. 



V. CONCLUSION 

In this article, we have characterized the potential energy surfaces (PESs) of the ozone 
molecule at large distances, where O and O2 weakly interact with each other, using two 
different approaches. Our calculations are based on the multipolar expansion of the elec- 
trostatic potential energy, which is expressed in terms of the electric properties of O and 
2 . 

Firstly, we have calculated the 27 asymptotic electronic PESs correlated to the 0( 3 P)- 
02(X 3 T,~) dissociation limit including spin-orbit couplings. All the PESs are found attrac- 
tive due to the dominant isotropic van der Waals interaction. For all spin-orbit states, the 
minimum energy is observed for a linear configuration of the three oxygen atoms. These 
long-range PESs can readily be connected at short distances to the most recent ones ob- 



tained by high-level ab initio calculations 



2l|, in order to derive global PESs for all values 



of the chosen internal coordinates. The overlap of the long-range and ab initio PESs on a 
large enough range of distances in the region of the LeRoy radius should avoid the tedious 
calculation of higher order terms in the multipolar expansion. 

Secondly, we have considered the interaction between an 2 (X 3 T l ~) molecule in a given 
spin-orbit state, vibrating in the ground level and rotating in a low level and an 0( 3 P) atom. 
We obtained one-dimensional potential energy curves (PECs) dominated by the isotropic 
van der Waals interaction. They are nearly parallel, which indicates that, when O and O2 
approach each other, the rotation of O2 is only weakly hindered by the presence of O at 
distances larger than the LeRoy radius. It is worthwhile to note that this hypothesis is 
generally assumed in models characterizingthe transition state, which is located at smaller 



distances, in the 0+0 2 exchange reaction [8 



55l | . However, such an hypothesis of weakly 
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hindered rotation is not always appropriate, as demonstrated for instance in the long-range 



45| 



interaction between a cesium atom and a ground state Cs2 molecule 

The one-dimensional PEC above can be smoothly matched to an arbitrary potential 
curve below the LeRoy radius so that vibrational energies and wavefunctions (along the 
O-O2 axis) can be calculated in the frame attached to the trimer (i.e. neglecting the mutual 
rotation of O and O2. They will be of relevance for estimating the density of vibrational 
levels of the O-O2 system close to the dissociation limit, namely with binding energy smaller 
than the remaining potential energy V(Rlr) ~ 30 cm -1 around the LeRoy radius. Such 
wavefunctions will be mostly concentrated at large distances, since the potential energy 
around the O3 global minimum is by far larger than V(Rlr). Once experimental energies 
are available near the dissociation limit, one could tune the potential energy in the inner 
zone in order to match the calculated energies to the observed one, and thus characterize 
the related energy levels with quantum numbers of the O-O2 system. 

Furthermore, such PECs could be of interest for the study of O-O2 collisions in the cold 
(T < 1 K) or ultracold (T < 1 mK) regimes. Several theoretical studies already aimed 



at modeling collisions involving oxygen 
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60], while several experimental investigations 



6l|-|63|. In 



suggested possible ways to create beams of cold atomic or molecular oxygen 
particular, we can expect such collisions to exhibit rotational unconventional isotopic effects. 
If the difference in rotational energy compensates the difference in vibrational zero-point 
energy, we obtain two quasi-degenerate channels of different isotopologues. In addition, the 
nuclear-spin symmetry can also play an important part in the rotational isotopic effects, 
as it forbids some rotational levels. For example, the channels 17 + 16 (^(A^ = 3) and 
16q + ie o 17 0(iV d = 4) on the one hand, and 18 + 17 18 0(N d = 3) and 17 + 18 2 (N d = 4) 
on the other hand, are only separated by 0.2 cm -1 ; but the second situation is not allowed, 
as 18 02 has no even rotational levels. 

It is worthwhile to mention that in contrast with similar studies devoted to low energy 
atom exchange reactions in atom-diatom systems (for instance AI+O2 64J, C+OH 65J, or 



O+OH 



33j), the adiabatic capture centrifugal theory [66[ conducted in the sudden approxi 



mation (ACCSA) 67| relying on the unique knowledge of the long-range asymptotic ground 
state potential energy of the atom-diatom system cannot be used in the ozone case, due to 
the presence of the transition state playing a role at distances lower than the LeRoy radius, 
and thus beyond the validity of the present approach. A calculation within this model would 



25 



undoubtedly lead to an overestimation of the atom exchange rate. 



ACKNOWLEDGMENTS 

We thank Pr. Vladimir Tyuterev for attracting our attention on the ozone problem. This 
work was done with the support of "Triangle de la Physique" under contract 2008-007T- 
QCCM (Quantum Control of Cold Molecules). 



[1] R. Schinke, S. Y. Grebenshchikov, M. V. Ivanov, and P. Fleurat-Lessard. Dynamical studies 
of the ozone isotope effect: A status report. Annu. Rev. Phys. Chem., 57:625-661, 2006. 

[2] M.H. Thiemens and J.E. Heidenreich. The mass-independent fractionation of oxygen: A novel 
isotope effect and its possible cosmochemical implications. Science, 219:1073, 1983. 

[3] K. Mauersberger. Ozone isotope measurements in the stratosphere. Geophys. Res. Lett., 
14:80-83, 1987. 

[4] S.M. Anderson, D. Hulsebusch, and K. Mauersberger. Surprising rate coefficients for four 

isotopic variants of 0+02+m. J. Chem. Phys., 107:5385-5392, 1997. 
[5] C. Janssen, J. Guenther, K. Mauersberger, and D. Krankowsky. Kinetic origin of the ozone 

isotope effect: a critical analysis of enrichments and rate coefficients. Phys. Chem. Chem. 

Phys. , 3(21):4718-4721, 2001. 
[6] K. Mauersberger, D. Krankowsky, C. Janssen, and R. Schinke. Assessment of the ozone isotope 

effect. Adv. At. Mol. Opt. Phys. , 50:1-54, 2005. 
[7] C. Janssen, J. Guenther, D. Krankowsky, and K. Mauersberger. Relative formation rates of 

50 o 3 and 52 o 3 in 16 o- 18 o mixtures. J. Chem. Phys. , 111:7179-7182, 1999. 
[8] Y.Q. Gao and R.A. Marcus. Strange and unconventional isotope effects in ozone formation. 

Science , 293:259-263, 2001. 
[9] Y.Q. Gao and R.A. Marcus. On the theory of the strange and unconventional isotopic effects 

in ozone formation. J. Chem. Phys., 116:137-154, 2002. 
[10] SY. Grebenshchikov, R. Schinke, P. Fleurat-Lessard, and M. Joyeux. van der Waals states 

in ozone and their influence on the threshold spectrum of O (x 1 ai). i. bound states. J. Chem. 

Phys. , 119:6512, 2003. 

26 



[11] D. Babikov. Entrance channel localized states in ozone: Possible application to helium nan- 
odroplet isolation spectroscopy. J. Chem. Phys. , 119:6554, 2003. 

[12] H.S. Lee and J.C. Light. Vibrational energy levels of ozone up to dissociation revisited. J. 
Chem. Phys. , 120:5859, 2004. 

[13] D. Babikov, B.K. Kendrick, R.B. Walker, R.T. Pack, P. Fleurat-Lesard, and R. Schinke. 
Metastable states of ozone calculated on an accurate potential energy surface. J. Chem. 
Phys. , 118:6298-6308, 2003. 

[14] S. Y. Grebenshchikov and R. Schinke. Towards quantum mechanical description of the un- 
conventional mass-dependent isotope effect in ozone: Resonance recombination in the strong 
collision approximation. J. Chem. Phy , 131:181103, 2009. 

[15] P. Fleurat-Lessard, S.Y. Grebenshchikov, R. Siebert, R. Schinke, and N. Halberstadt. Theo- 
retical investigation of the temperature dependence of the 0+02 exchange reaction. J. Chem. 
Phys. , 118:610, 2003. 

[16] M.V. Ivanov and D. Babikov. Efficient quantum-classical method for computing thermal rate 

constant of recombination: Application to ozone formation. J. Chem. Phys., 136:184304, 2012. 
[17] R. Siebert, R. Schinke, and M. Bittererova. Spectroscopy of ozone at the dissociation threshold: 

Quantum calculations of bound and resonance states on a new global potential energy surface. 

Phys. Chem. Chem. Phys. , 3:1795-1798, 2001. 
[18] R. Siebert, P. Fleurat-Lessard, R. Schinke, M. Bittererova, and S.C. Farantos. The vibrational 

energies of ozone up to the dissociation threshold: Dynamics calculations on an accurate 

potential energy surface. J. Chem. Phys. , 116:9749-9767, 2002. 
[19] P. Rosmus, P. Palmieri, and R. Schinke. The asymptotic region of the potential energy surfaces 

relevant for the o (p)+o 2 {x 3 (Tg )^o 3 reaction. J. Chem. Phys. , 117:4871-4877, 2002. 
[20] F. Holka, P.G. Szalay, T. Mtiller, and V.G. Tyuterev. Toward an improved ground state 

potential energy surface of ozone. J. Phys. Chem. A, 114:9927-9935, 2010. 
[21] R. Dawes, P. Lolur, J. Ma, and H. Guo. Communication: Highly accurate ozone formation 

potential and implications for kinetics. J. Chem. Phys. , 135:081102, 2011. 
[22] M. Tashiro and R. Schinke. The effect of spin-orbit coupling in complex forming o( 3 p)+02 

collisions. J. Chem. Phys. , 119:10186, 2003. 
[23] B. Bussery-Honvault, F. Dayou, and A. Zanchet. Long-range multipolar potentials of the 18 

spin-orbit states arising from the C ( 3 P) + OH(X 2 n) interaction. J. Chem. Phys. , 129:234302, 



[24 
[25 
[26 
[27 
[28 
[29 
[30 



[31 

[32 
[33 



2008. 

B. Bussery-Honvault and F. Dayou. Si ( 3 -P) + OH (X 2 II) interaction: Long-range multipolar 

potentials of the eighteen spin-orbit states. J. Phys. Chem. A, 113(52):14961, 2009. 

M. Lepers, O. Dulieu, and V. Kokoouline. Photoassociation of a cold atom-molecule pair: 

long-range quadrupole-quadrupole interactions. Phys. Rev. A, 82:042711, 2010. 

R.J. LeRoy. Long-range potential coefficients from RKR turning points: Cq and Cg for 

J B( 3 n Gu +)-state Cl 2 , Br 2 , and I 2 . Can. J. Phys. , 52:246, 1974. 

M. Bartolomei, E. Carmona-Novillo, J. Campos-Martinez M.I. Hernandez, and R. Hernandez- 
Lamoneda. J. Comput. Chem. , 32:279-290, 2010. 

M. Tinkham and M. W. P. Strandberg. Theory of the fine structure of the molecular oxygen 
ground state. Phys. Rev. , 97:937-951, 1955. 

D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii. Quantum theory of angular 
momentum. Leningrad, 1988. 

D. Spelsberg, T. Lorenz, and W. Meyer. Dynamic multipole polarizabilities and long range 
interaction coefficients for the systems H, Li, Na, K, He, H~, H 2 , Li 2 , Na 2 , and K 2 . J. Chem. 
Phys. , 99:7845, 1993. 

P.W. Langhoff and M. Karplus. Pade approximants for two-and three-body dipole dispersion 
interactions. J. Chem. Phys. , 53:233-250, 1970. 
EPAPS no. NNN. 



T. Stoecklin, B. Bussery-Honvault, P. Honvault, and F. Dayou. Asymptotic potentials and 
rate constants in the adiabatic capture centrifugal sudden approximation for X + OH(A 2 n) — > 
OX + R( 2 S) reactions where X = 0( 3 P), S( 3 P) or N( 4 5). Comp. Theor. Chem. , 990:39-46, 
2012. 

[34] G.L. Gutsev, P. Jena, and R.J. Bartlett. Chem. Phys. Lett. , 291:547, 1998. 

[35] A.K. Das and A.J. Thakkar. J. Phys. B , 31:2215, 1998. 

[36] M. Medved, P.W. Fowler, and J.M. Hutson. Mol. Phys. , 98:453, 2000. 

[37] D. B. Lawson and J. F. Harrison. Distance dependence and spatial distribution of the molec- 
ular quadrupole moments of H 2 , N 2 , 2 , and F 2 . J. Phys. Chem. A , 101(26):4781-4792, 
1997. 

[38] A. Kumar, W. J. Meath, P. Biindgen, and A. J. Thakkar. Reliable anisotropic dipole prop- 
erties, and dispersion energy coefficients, for o evaluated using constrained dipole oscillator 

28 



strength techniques. J. Chem. Phys. , 105:4927-4937, 1996. 
[39] D.A. Dummur A.D. Buckingham, R.L. Disch. J. Am. Chem. Soc , 90:3104, 1968. 
[40] M.A. Buldakov, I.I. Ippolitov, B.V. Korolev, I.I. Matrosov, A.E. Cheglokov, V.N. 

Chrerepanov, Yu.S. Makushkin, and O.N. Ulenikov. Spectrochim. Acta A, 52:995, 1996. 
[41] E. Cohen and Birnbaum G.J. J. Chem. Phys. , 66:2443, 1977. 
[42] N.J. Bridge and A.D. Buckingham. Proc. R Soc. London Ser. A , 295:334, 1966. 
[43] M. Lepers, R. Vexiau, N. Bouloufa, O. Dulieu, and V. Kokoouline. Photoassociation of a cold 

atom-molecule pair: second-order perturbation approach. Phys. Rev. A, 83:042707, 2011. 
[44] M. Lepers and O. Dulieu. Ultracold atom-dimer long-range interactions beyond the l/r n 

expansion. Eur. Phys. J. D , 165:113-123, 2011. 
[45] M. Lepers and O. Dulieu. Long-range interactions between ultracold atoms and molecules 

including atomic spin-orbit. Phys. Chem. Chem. Phys., 13:19106-19113, 2011. 
[46] A. Van der Avoird and G. Brocks. The O2-O2 dimer: magnetic coupling and spectrum. J. 

Chem. Phys , 87:5346-5360, 1987. 
[47] M.-L. Dubernet and J. M. Hutson. Atom-molecule van der Waals complexes containing open- 
shell atoms, i. general theory and bending levels. J. Chem. Phys. , 101:1939-1958, 1994. 
[48] T. Y. Chang. Moderately long-range interatomic forces. Rev. Mod. Phys., 39(4):911-942, 

1967. 

[49] D. M. Bishop. Molecular vibrational and rotational motion in static and dynamic electric 

fields. Rev. Mod. Phys. , 62(2):343-374, 1990. 
[50] A.C. Allison, S.L. Guberman, and A. Dalgarno. A model of the Schumann-Runge continuum 

of 2 . J. Geophys. Res. , 91(A9):10193-10, 1986. 
[51] R.S. Friedman. Oscillator strengths of the schumann-runge bands of isotopic oxygen molecules. 

J. Quant. Spec. Rad. Trans. , 43:225-238, 1990. 
[52] J. Wang, D.G. McCoy, A.J. Blake, and L. Torop. Effects of the close approach of potential 

curves in photoabsorption by diatomic molecules-ii. temperature dependence of the O2 cross 

section in the region 130-160 nm. J. Quant. Spec. Rad. Trans. , 38:19-27, 1987. 
[53] R.B. Singh and D.K. Rai. Potential-energy curves for 0^+, N^, and CO + . J. Mol. Spec , 

19:424-434, 1966. 

[54] D. Spelsberg and W. Meyer. Ab initio dynamic dipole polarizabilities for O2, its photoab- 
sorption spectrum in the schumann-runge region, and long-range interaction coefficients for 

29 



its dimer. J. Chem. Phys. , 109:9802-9810, 1998. 
[55] B.C. Hathorn and R.A. Marcus. An intramolecular theory of the mass-independent isotope 

effect for ozone, ii. numerical implementation at low pressures using a loose transition state. 

J. Chem. Phys. , 113:9497-9509, 2000. 
[56] A. V. Avdeenkov and J. L. Bohn. Ultracold collisions of oxygen molecules. Phys. Rev. A, 

64:052703, 2001. 

[57] A. Volpi and J.L. Bohn. Magnetic-field effects in ultracold molecular collisions. Phys. Rev. 
A, 65:052712, 2002. 

[58] K. Tilford, M. Hoster, P.M. Florian, and R.C. Forrey. Cold collisions involving rotationally 
hot oxygen molecules. Phys. Rev. A , 69:052705, 2004. 

[59] J. Perez-Rios, M. Bartolomei, J. Campos-Martmez, M.I. Hernandez, and R. Hernandez- 
Lamoneda. Quantum- mechanical study of the collision dynamics of 02( 3 S~) + 02( 3 S~) on a 
new ab initio potential energy surface. J. Phys. Chem. A, 113:14952, 2009. 

[60] J. Perez-Rfos, J. Campos-Martmez, and M.I. Hernandez. Ultracold O2 + O2 collisions in a 
magnetic field: On the role of the potential energy surface. J. Chem. Phys., 134:124310, 2011. 

[61] D. Patterson and J.M. Doyle. Bright, guided molecular beam with hydrodynamic enhance- 
ment. J. Chem. Phys. , 126:154307, 2007. 

[62] E. Narevicius, A. Libson, C.G. Parthey, I. Chavez, J. Narevicius, U. Even, and M.G. Raizen. 
Stopping supersonic oxygen with a series of pulsed electromagnetic coils: A molecular coilgun. 
Phys. Rev. A , 77:051401, 2008. 

[63] O. Bucicov, M. Nowak, S. Jung, G. Meijer, E. Tiemann, and C. Lisdat. Cold SO2 molecules 
by Stark deceleration. Eur. Phys. J. D , 46:463, 2008. 

[64] D. Reignier, T. Stoecklin, S.D. Le Picard, A. Canosa, and B.R. Rowe. Rate constant calcu- 
lations for atom-diatom reaction involving an open-shell atom and a molecule in a a elec- 
tronic state application to the reaction Al( 2 P 1 /2 >3 /2)+02(X 3 I!-) -)• A10(X 2 I! + )+0( 3 .P2,i,o)- 
J. Chem. Soc, Faraday Trans. , 94(12):1681-1686, 1998. 

[65] M. Jorfi, B. Bussery-Honvault, P. Honvault, T. Stoecklin, P. Larregaray, and P. Halvick. 
Theoretical sensitivity of the C( 3 P)+OH(A 2 n)^CO(A 1 S+)+H( 2 S') rate constant: The role 
of the long-range potential. J. Phys. Chem. A , 114(28):7494-7499, 2010. 

[66] R.D. Lcvinc and R.B. Bernstein. Molecular Reaction Dynamics and Chemical Reactivity . 
Oxford University Press, 1987. 



30 



[67] D.C. Clary. Rates of chemical reactions dominated by long-range intermolecular forces. Mol. 
Phys. , 53(1):3-21, 1984. 



31 



